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1 Introduction 



Motivated by the experimental progress in relativistic heavy ion colhsions and the 
recreation of the quark gluon plasma, several questions have received substantial at- 
tention in the past few years. What happens to hadrons in the deconfined quark-gluon 
plasma? Do bound states persist? What is rate of photon and dilepton production 
from a hot QGP? How effectively are energy-momentum and charge transported? 
How long, or rather how short, are the typical relaxation times for hydrodynamic 
fluctuations? 

Since this information is encoded in spectral functions, it is prohibitively difficult 
to access it directly from euclidean correlators obtained with lattice QCD, due to 
the intricacy of performing the analytical continuation from imaginary to real time. 
However, recent progress has been made by applying the Maximal Entropy Method 
(MEM) jTj to this problem. An (incomplete) list of high temperature studies includes 
the possible survival of hadronic bound states in the deconfined quark-gluon plasma 
I21IS1I1]; thermal dilepton rates and transport coefficients [Oj^ 

In a spectral function investigation, the low-energy region < T is of particular 
interest, since it is expected to be the most affected by nonperturbative medium ef- 
fects. However, the reconstruction of spectral functions at small energies <C T is 
hindered by the insensitivity of euclidean correlators to details of spectral functions 
at these energies [S]. This is especially important for the calculation of transport 
coefficients where by definition the interest is in the limiting value of current-current 
spectral densities as — 0. Experience with the reconstruction of spectral densities 
in the low-energy region can be obtained by studying the simpler (but still nontrivial) 
problem of meson spectral functions at nonzero momentum above the deconfinement 
transition. Due to e.g. the scattering of quarks with gauge bosons below the light- 
cone (Landau damping), these spectral functions are expected to have a nontrivial 
structure. Since in the confined phase one expects to find mesons moving relative 
to the heatbath, described by simple quasiparticle spectral functions, increasing the 
temperature from below to above the transition temperature should result in a drastic 
change in those spectral functions. 

Our aim in this paper is to provide a reference point for such an analysis on the 
lattice in the infinite temperature limit. It is therefore similar in spirit as Ref. jTU], in 
which a study at zero momentum was performed. The paper is organized as follows. In 
the next section we give complete analytical expressions for continuum meson spectral 

^We note here that Ref. |7j does not use an MEM analysis, but instead employs an Ansatz which 
was proposed in Ref. [H] and criticized in Ref. [H]. 
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functions at nonzero momentum and fermion mass in the infinite temperature limit 
and discuss several features. In Section |21 we derive simple expressions for meson 
spectral functions for Wilson and staggered lattice fermions. We briefly comment on 
the value of the euclidean correlator at the midpoint and on the use of local and point 
split currents. The main results are shown in Section |3J where we contrast spectral 
functions obtained with Wilson and staggered fermions with the continuum results. 
Section El contains a short summary. 



2 Continuum 

We consider meson spectral functions with quantum numbers if, defined as 

p^,(t,x) = ([J^,(t,x), 4(0,0)]), (1) 

with Jh{t,x.) = g(r, x)r//g(r, x) and Th = {1, 75, 7^^, 7'^75}.^ They are related to 
euclidean correlation functions, 

GH(r,x) = (JH(r, x)4(0,0)), (2) 

via the standard integral relation 

Gnir, P) = ^ ^ ^'(^' ^)Ph{^, p), (3) 

with the kernel 

where nsico) = l/^e^^"^ — 1) is the Bose distribution. At lowest order in the loop 
expansion, the euclidean correlators read in momentum space^ 



Gh{p) = - ^^ti s{K)rHS{p + K)-f'r^Hi' 



(5) 



where P = {iujn, p) with Un = 27mT {n G Z) the Matsubara frequency in the 
imaginary-time formalism, and 



K 



^In this section the gamma matrices obey 7"^ ~ 7*^, 7't =^ —7*, and 7^ = 75. The anticommuta- 
tion relations are {7^,7*^} = 25'^" and {7'^, 75} = with g^^'^ = diag(+, — , — , — ). 
'^7''r|^7*' appears since the original correlator is of the form {JJ^), not {J J). 
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The fermion propagators are given by 



S{K) = -\ = - r ^P^^, (7) 



iUn'y^ — 7 ■ k — m 27r zu;„ — u ' 

where a;„ = {2n + l)7rT (n G Z) is a fermionic Matsubara frequency and Pf{^i k) the 
spectral density of the fermion, 

Pf{K) = + m) p{K) = + m) 2ti sgn(A;°)5(A;2 - ^2)^ (g) 



with ciJk = vkM-?7?. 

Using the spectral representation for the fermion propagators, it is straightforward 
to arrive at 

Ph{P) = 2lmGH{iuJn ^uj + ?0+, p) 

= N, f tr + m)TH{j^ + m)7°rl,7° p{K)p{R) [ripik') - np(r°)] , (9) 

with P = (a;, p), R = P + K and npioj) = l/(e'^/^ + 1) is the Fermi distribution. 

To facilitate the comparison with the lattice expressions below, we give here the 
result with the integral performed. 



k ■ r m 



(1) ^(2)K-I (3) 
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a\j I [1 — nF{uJ\^ — nF{uJr)\ 5{uj — CJk — ^^r) 



-{uj^-uj)y (10) 

The first line corresponds to scattering and contributes only below the lightcone (a;^ < 
p^. Landau damping), while the second line corresponds to decay, contributing above 
threshold (u;^ > + 4m^). The coefficients af^ arise from the three nonzero traces 
over the gamma matrices in Eq. (jH)) and depend on the channel under consideration. 
They are listed in Table ^ 

The remaining integrals can be performed as well. In terms of 



1 / 4to2 
P± = -[uj±pP{P)], PiP) = \l , s = u;'-p^ (11) 
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i H 




a^j 




Ps 


1 


1 


-1 


1 


Pps 


75 


1 


-1 


-1 






1 


1 


1 


pii 


7^ 


3 


-1 


-3 


Pv 




2 


-2 


-4 


pT 


7°75 


1 


1 


-1 




7*75 


3 


-1 


3 


PA 


7^75 


2 


-2 
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Table 1: Coefficients afj for free spectral functions in different channels H. In the 
case of 7* and 7*75, the sum is taken over z = 1, 2, 3. We defined py = —g^vP^^ and 
pA = -g,,upt''- 



the final expression in the continuum reads 



Ph{P) = 0(s - 4m^ 



NrT"^ 



TT 



24T2 
1 



(1) , ^^2 , ,2^2^D^^ ^(2) _ ^_2^(3)1 . 1+6"^'+/^ 

/T 
P-/T 



ApT 



2T 



In 



1 + e~P-/T 



In 



TT 

1 + e-P+/^ 
1 + eP-/^ 



+ — [Li3(-e-^+/n-Li3( 
p 



{u' - p'p\P)) a^£> + {p' - u'P^P)) a^^ - 4mV^ 



ApT 



(2) 



.2^(3) 



+ {o!^S + (^^^) [Li2(-e-^>/^) - Li2(-e^--/^ 

+— [Li3(-e-^"+/^)-Li3(-e^^-/^ 
p 



■ (12) 



We now discuss several features. First consider the asymptotic behaviour at large uj. 
We find that all spectral functions increase with cj^, 



lim Ph{.P) = ©(-s - 



,2 



in 



(13) 



as expected from naive dimensional arguments, except when there is a cancellation. 
This happens for Th = 7'',7''75, for which we find instead 



lim p"'^(P) = e(s - 4m')—/, 
lim pf(P) = e(s - Am^)— ip" + 6m2) 



(14) 



For the vector current this behaviour can be understood from current conservation 
d^j^ = 0. Since at large uj the effect of finite temperature is exponentially suppressed, 
we may use the zero temperature decomposition. 



/7^^(P) = 2Imn^^(P) = 2 (PV - P^P"") Imn^(P^ 



(15) 



which explains the behaviour above. Current conservation also relates the other 
components of p^^ , 

ujp^^P) =p'p'\P), ujp^'iP) =p'p'^{P), cuV°(P) =pVp'^(P). (16) 

Since the axial vector current is not conserved, 

d^jt^ = 2mj5 + anomaly, (17) 

similar relations do not hold for pg''. However, in the free case considered here, we 
find 

u;V5°(^,0) = 4mVs(^,0). (18) 

Any deviation from this is therefore due to the U(1)a anomaly. 
In the zero momentum limit, the spectral functions reduce to"^ 



N, 



Ph{uj, 0) = e(cj2 - 4m2)^Vcu2 -4m2 [1 - 2nF{uj/2)] 



iTTLU 



UJ a 



+Am a 



with 



''H 



h 



2'ituj6{uj)Nc 



''H 



''H 



''H 



(19) 



(20) 



a^"* — a'"^' = 2bH, except for the axial currents {AQ,Ai), where we find bn = (1,-2) instead of 
bn = (0,3). Note that the coefficients in Ref. ^01 disagree with relation H18|l . Note also that the 
normalization differs by a factor of 27r and that the overall signs for ps and p.s are opposite. 



*A comparison with the coefficients an and bn in Table 2.1 of Ref. ^U] yields a^-* — a^' = 2aH 



,(3) 
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In the massless case Ii = I2 = T'^/Q. The term proportional to ujS{uj) is all that 
remains from the scattering contribution below the lightcone in Eq. It gives a r 

independent contribution to the euclidean correlator since the kernel K{t,lj) ~ 2T/uj 
for small u. In particular, charge conservation dictates the form of and at 
zero momentum, 

0) = 27rxuj6iu;), G°°(r, 0) = Tx, (21) 

which is not altered by interactions, although the value of the charge susceptibility 
X is. At the order computed here, x = '^NJi- For completeness we give here the 
euclidean correlator at zero momentum and mass 



jv 



6 



(1) (2) , 3 / (1) _ (2)\ Su + u cos(2m) - 2 sin(2M) 
/V / sm u 



(22) 



where u = 2-kT{t - 1/2T). 

Finally, it follows from the spectral decomposition 

Ph{P) = \Y.\ WHmm)\\2nY5\P + P„ - P„) (e"^"/^ - e^^™/^) , (23) 

n,m 

where Z is the partition function, that all spectral functions for a single current Jh are 
odd and positive semi-definite for positive argument, i.e. upnioj^Tp) > 0. Obviously, 
spectral functions that are defined as the difference between such spectral functions, 
such as pv = ~ and pA = P5 — can turn negative. Indeed, it is easy to 
see that pv('-^5P) is negative for small uo ii < 2m?. All other spectral functions 
increase linearly with uj for small u and nonzero p. 

Although not the topic of this paper, we briefly mention how corrections due 
to interactions appear at very high temperature. First of all, for soft momentum 
IpI ~ gT, a hard thermal loop [TTj calculation is needed, see e.g. Refs. |T^ITHlin] 
for such studies. The gap in the spectrum for < uj"^ < + Am? is filled when two 
loop diagrams are included, due to e.g. bremsstrahlung jT^j. Around the lightcone 
the loop expansion breaks down due to the Landau-Pomeranchuk-Migdal effect and 
an infinite series of ladder diagrams contribute at leading order in the strong coupling 
constant P^ llTj. Finally for very soft momenta and energies, the structure of current- 
current spectral functions is determined by general hydrodynamical considerations 
|18j . So far a diagrammatic calculation in this regime has been carried out only 
in the case of the spatial vector spectral function p**(a;,0) in the limit of exactly 
zero momentum and vanishing energy a; ^ 0, which is relevant for the electrical 
conductivity: see Refs. [121 120] for details on the weak coupling result at leading- 
logarithmic order and Ref. [25 for the large Nf result. 
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3 Lattice 



3.1 Wilson fermions 



In this section we derive expressions for meson spectral functions on a lattice with 
A'^ X N-r sites. The lattice spacing is denoted with a in the spatial directions and with 
ttr in the temporal direction, ^ = a/ar is the anisotropy parameter. The temperature 
is related to the extent in the imaginary time direction, T = l/(A^T-aT-). We start 
with standard Wilson fermions. The lattice fermion propagator (with coefficients 
r^pace = r) reads^ 



S{K) 



—274 sin k/i — z/Ck + r4 (1 — cos k^) + M.^ 
sin^ /C4 + /C^ + [r4 (1 - cos k^) + M\,f 



where 



1 

/Ck = - ^ 7j sin ki 



(1 — cos kj) + m 



(24) 



(25) 



We use periodic boundary conditions in space, ki = 2Tmi/N„ with Ui = —Nfj/2 + 
1, —Nrj/2 + 2, . . . , Ncr/2 — 1, Ncr/2 for i = 1,2,3, and antiperiodic boundary conditions 
in imaginary time, ^4 = 7r(2n4 + 1) /N^ with m = -Nr/2 + 1, -Nr/2 + 2, . . . , Nr/2 - 
l,iV,/2. 

To make a smooth connection with the expressions in the continuum we follow 
Ref. PU] and use the mixed representation of Carpenter and Baillie 



5(r, k) = 74^4(r, k) + J] 7.5.(r, k) + ]15„(r, k). 



(26) 



In order to avoid the doubler in the time direction, we proceed with r4 = 1, so that 
(for < r < l/T) 



5*4 (r, k) = 5*4 (k) cosh(r-E'k) 
5,(r,k) = S',(k) sinh(fEk). 

Suij, k) = S'„(k) sinh(rEk) 



2(1 + Mi 



(27) 



^In this section the gamma matrices are hermitian, 7^, =7^, 75 =75, and obey {7^^,71/} = 2(5^i,, 
{7^, 75} = 0. They are related to the gamma matrices of the previous section as 74 = 7*^, 7^ = —17'. 
We use lattice units a = 1. 
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Here r = r — 1 / 2T and 



sinh 

28], cosh(Ek/2T) ' 
1 i sin /cj 
^2£:kCosh(Ek/2r)^ 



c. ^ _ 1-cosh {E^/O + M^ 
2£:kCosh(Ek/2r) ' 



with £'k = (1 + A^k) sinh {E\,/^). The single particle energy E'k is determined by^ 

The final term in 5'„(r, k) is the sole remnant of the nonpropagating time doubler; 
below we consider < r < 1/T. The propagator satisfies S{—t, k) = '~^^S\t, k)75. 
The correlators we are interested in are of the form 

Gh{t, p) = 5^ tr S{r, \<i)THS{-T, v)Th, (30) 

k 

where again r = p + k. Inserting Eq. (|2b|l gives the euclidean correlator^ 
Gh{t. p) = ^ E [«h ^4(r, k)5l(r, r) - ag^ ^ S,{t, k)4(r, r) 

k i 

-ag)^„(r,k)5t(r,r)], (31) 

where the coefficients o'^ are the same as before (see Table CJ. 

We will now extract the spectral functions in a form that closely resembles the con- 
tinuum expressions. In the terms 55^ we encounter products of hyperbolic functions. 
These can be written as 

sinh(r-E'k) sinh(r£'r) = 7 / duo cosh(co'r) 6{uj — E\, — Ej.) + S{uj + -Ek + Ej.) 

^ J —00 

-6{UJ - Ek + ^r) - S{UJ + Ek - ^r)l ,(32) 



^The factor 1/^ is included so that in the continuum hmit Ey^ ^ Wk = V^ + m^ (with a = 1). 
^Note that we now start from {Jh{t,x.)Jh{0,O)) rather than from Eq. 0. This only affects the 
overall sign in some channels, which has been adjusted to agree with the continuum one. 
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and similarly for the product of two hyperbolic cosines. Noting that the factor 
cosh(ci;'f) is the sole place with r dependence and that it is of the same form as 
in the kernel (jlj), it is straightforward to write the above expression for Gh{t, p) as 



Jo 

and read off the expressions for the lattice spectral functions, 



(33) 



2TJ 



a^i^S,{k)Sl{r) + ag) ^ 5.(k)5l(r) + ag)^„(k)5t (r 

i 

a«54(k)^l(r) - ag) 5^^,(k)^t(r) _ ag)5„(k)5t(r) 



+ 



6{uo -E^- E^) 



(34) 



This result can be directly compared with the continuum expression using Eq. 
(j2Hl) and realizing that 



sinh(cu/2T) 



2 [npiE^) - nF{E,)] 



if LJ = Er-E- 



k, 



cosh(Ek/2r) cosh(Er/2T) [2 [1 - npiEi,) - upiEr)] ii u = E^ + E^. 
3.2 Staggered fermions 



(35) 



In the case of staggered fermions we perform the analysis with naive fermions, since 
this leads to equivalent results [22] • Taking therefore r4 = rspace = yields the fermion 
propagator with 

^4(r,k) = S'4(k) (1 - (-1)-/"^) cosh(f£;k), 
^,(r,k) = ^.(k) (1 + i-iy/^^) sinh(f£;k), 

5„(r, k) = SM (1 + (-1)-/'^^) sinh(f Ek), (36) 

where S'4(k) and S'j(k) are as in Eq. (PHj) with = cosh (-Ek/O ^inh (i^k/O; 

1 



Su{t, k) 



'^2^:kCosh(£;k/2T)' 
10 



(37) 



The single particle energy Ey_ is now determined by 



cosh {EJO = y/l + /C2 + (m/02. (38) 

Using the same steps as before, the euchdean meson correlator takes again the form 
(j3ip and can be written in a spectral representation as 

G^(r, P) = 2 ^K{t, u) [p-'™(a;, p) - {-ly/^^ p^^^u, p)] , (39) 

with the same kernel as above. 

The desired spectral function p^'^°(u;,p) is exactly as in Eq. (jSl|), whereas the 
staggered partner p^'™(ci;,p) has the same form but with coefficients a^"* = 

= — a^H — ~^h^ ■ This staggered contribution pH represents the spectral 
function in the channel related to the original pn by replacing Th = 'j^'^dTn 

|24j . Note that in particular the pseudoscalar (scalar) spectral function mixes with 
the zero'th component of vector (axial vector) current spectral function. In an actual 
MEM investigation, the staggered partners can be disentangled using an independent 
analysis on even/odd timeslices, which yields the linear combinations p^'™(u;,p) =i= 
-naivej^^ Finally, in order to compare the naive lattice spectral functions with the 
continuum and the Wilson ones, we divide pjf"'^ by a factor of 8, which takes care of 
the space doublers. 

3.3 Midpoint of the euclidean correlator 

In the midpoint r = 1/2T (f = 0), the hyperbolic functions in the fermion propagator 
S{t, p) take simple values, and it is easy to see that 

Gh{1/2T,p) = ^^aiJ)54(k)5l(r). (40) 

k 

This implies that the channel dependence of the value at the midpoint enters only 
via a^^\ 

Analogous expressions hold for naive fermions and in the continuum, so that one 
can write 

G^(l/2T,p) = ai;)C(p), (41) 
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with 



C'continuum (p) -^c I r"7 

ik cosh(u;i, 



1 



3h(cUk/2T) cosh(cJr/2T) 

r ri = ^v - - 

WiisonlPj ^3 (1 + cosh(£;k/2T) (1 + Mr) cosh{Ej2T) ' 

odd/ X _ 4iVc 1 1 

•-"naive IPJ — ~7F 



L3 ^ cosh(Ek/0 cosh(Ek/2T) cosh(Er/0 cosh(Er/2r) ' 



^^n'll'^™^(p)=0. (42) 

Combining this resuh with the relation between the euchdean correlator and the 
spectral function pn in Eq. Q, yields a constraint for the free spectral density 

««(V2r,p)^f ^-^^a-C(p). (43) 

In the case of naive fermions this gives 

G-'™(1/2T, p) = 2 r ^ . ,/ [p?f^^<=(a., p) T pT™(^, p)] , (44) 
Jq 2n sinh(co'/2r) -' 

for Nr/2 even/odd, from which we find 

/•°^ du pgf-(u;,p) ^ /-"^ rfo; pTI^ ^ 1(1)^^^/2 odd. . .... 

Jo 27rsinh(cj/2T) 27r sinh(cu/2T) 4"^^ ""^-ve IPJ- I 

Although the free spectral functions in the various channels are distinctly different, 
we conclude that the integral of ph{uj,p)/ sinh(a;/2T) is in all cases related to C(p) 
given above, both in the continuum and on the lattice. 

3.4 Point split current 

The local vector current we have considered so far, = ip'Juip, is not exactly conserved 
on the lattice. Instead, the conserved current is 

3^{x) = tpix + a/,)P+V^(x) - tlj{x)p-tlj{x + a^), (46) 

with = |(r^ ± 7^), Ti = r and d^ = fta^. Here we present a short analysis 
comparing the two currents. 
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A correlator especially sensitive to the difference between the local and the con- 
served current is G^^{t,0) at vanishing momentum, since its form is determined by 
charge conservation, see Eq. (PT|) . In particular it should be r independent. On the 
lattice, the correlator for the local current is 



G°°(r, 0) = ^J2[ - E - \Suir, k)f]. (47) 



k 

After some algebra, this can be written as 

L3 4^ (1 



k 



A^k)'cosh'(Ek/2r) 



0)-^^-V l-(-ir/--cosh(2fEk) . . 

^naivel^.^J ^3 cOsh2(Ek/e) COsh2(Ek/2T) ' ^ ^ 

For Wilson fermions we find a r independent result at leading order. However, since 
the local current is not related to a symmetry, dependence on r is expected to arise 
when interactions are present. This is easy to study in actual simulations. For naive 
fermions we indeed find a r dependent result. 

With the conserved current the situation should be different. We find for Wilson 
fermions, with = 1, 

<^wiison(^> 0) = {S^ir - ctr, k)S'l(r + a^, k) - S'„(r - a^, k)Sl{T + a^, k) 



k 



J2 \S^ir, k)p - S,iT - a,, k)SliT + a., k) + 5„(r - a., k)^l(r + a., k)}, (49) 



and for naive fermions 
2N, 



GZUr, 0) = jfYl {^4(r - ar, k)Sl{r + a., k) + \S,{t, k)\' - |5„(r, k)!^ 

k 

-Suir - a^, k)SiiT + a., k) - [5,(r - a., k)^/(r + a., k) + |S',(r, k)!^] }. (50) 

i 

Indeed, this yields the anticipated result for a conserved current. 
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In both cases G'°°(r, 0) is now r independent; this should remain to be the case when 
interactions are included. Moreover, the lattice susceptibility takes the same form 
as in the continuum, see below Eq. (PT|) . The factor 1/2 in the naive case appears 
because of the contribution from the time doublers. 

If one is interested in the reconstruction of vector spectral functions for e.g. thermal 
dilepton production jS], it may be important to use the properly conserved current. 
It would therefore be interesting to compare spectral functions obtained with local 
and point split currents in the interacting case. 

4 Comparison 

We now contrast the meson spectral functions obtained for free Wilson and staggered 
lattice fermions with the continuum ones. The lattice meson spectral functions ob- 
tained above can be analysed for finite A'^,- and A'^. For small the discreteness 
inherent in the definition of spectral functions (see e.g. the spectral decomposition 
(|23p ) is clearly visible. Following Ref. jTUj we therefore take the thermodynamic limit 
N„ oo and focus on the effect of finite A'^.^ In all figures the nonzero external 
momentum p = 4T and the fermion mass m = T. For Wilson fermions we show 
results with r = 1. The anisotropy parameter ,^ = 1, except in the bottom part of 
Fig. El We only show meson spectral functions obtained with local operators. 

In Fig. n we show the scalar and pseudoscalar spectral functions for Wilson (left) 
and staggered (right) fermions. In order to emphasize the effects of the lattice cutoff, 
Pps^s is divided by u'^ in the top figures. The continuum result then reaches a constant 
value (3/4-7r) for large u, see Eq. p3|) . Instead, on the lattice there is a maximal 
energy cUmax, determined by the delta function 6{u! — E]^ — Ey.). Since the external 
momentum p is small with respect to momenta at the edge of the Brillouin zone, the 
maximum value for Wilson fermions (with r = 1) is determined by fermion momenta 
k = (tt/o, vr/a, vr/a) [TP], which gives 



and for staggered fermions by fermion momenta k = (7r/2a, 7r/2a, 7r/2a), which yields 

*Iii practice we take ^ 1500—2000, replace the delta functions in Eq. H34|l with block functions 
with width Aw and height l/Aw, and divide the lo interval in bins. We used = 1000. The 
bin width is determined by Aw = uiraax/N^ where Wmax is discussed below. See also |10) . 




(52) 
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Figure 1: Pseudoscalar and scalar spectral functions pps,s(i^, p)/t^^ (above) and 
Pps,s(i^, p)/7'^ (below) as a function of uj jT for A''^ = 16, p/T = A,m/T = 1 and 



.naive 
max 



2iV.ln Ve + i + a'm^ + VS + a'm^^ ^^^^ 

The maximum value is smaller for staggered fermions. The cusps in the plots originate 
from the corners of the Brillouin zone. Both for continuum and staggered fermions, 
we find that the scalar and pseudoscalar channel are indistinguishable for large u. 
The reason is that the finite fermion mass is negligible for such large energies. In the 
case of Wilson fermions the Wilson mass term breaks the chiral symmetry completely 
and the scalar and pseudoscalar spectral functions differ. 

The spectral functions vanish for energies p < uj < \/p^ + Am?. The physically 
interesting contribution below the lightcone appears as a divergent one in the top 
plots. We therefore show pps,s/^^ in the plots on the bottom. The spectral functions 
increase linearly for small uj and vanish at the lightcone. Due the finite fermion 
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Figure 2: Effect of finite lattice spacing on the 'lightcone' and 'threshold' for 
fixed pL ^ 25 and m = p/2. The continuum values are u/p = 1 resp. u/p = 
^1 + 4:m^/p'^ = V2. 



mass, the scalar and the pseudoscalar channel are now physically distinct. The main 
lattice artefact in this region appears to be the mismatch between the location of 
the lightcone in the continuum and the lattice theory. This is due to the difference 
between continuum and lattice dispersion relations. To study this further, we define 
the lattice 'lightcone' and 'threshold' via 

lightcone: = max(£'k — £'k+p), threshold: u; = min(£'k + -E'k+p)- (54) 

k k 

In Fig. 121 we show the result as a function of 1/N„ ~ a for fixed momentum p = 
(0, 0, Svr/aiVo-) {pL f» 25) and m = p/2. As expected, the continuum and lattice 
results agree for decreasing l/N^j (decreasing lattice spacing), but for finite a the 
corrections can be substantial. 

The effect of increasing is demonstrated in Fig. El for the pseudoscalar spectral 
function pps for fixed ^ = 1 (top) and for the vector spectral function p** for fixed 
^/Nt- = aT (bottom). As expected from Eqs. and ci^max increases with 

A^^. In the anisotropic large A^,- seems to lead to a better improvement for 

Wilson than for staggered fermions. In Fig. |31 we present our results for As we 
emphasized in Eq. (fT^ . due to current conservation this spectral function does not 
increase with a;^ for large u, but instead reaches a constant value p'^/2'K. This can 
indeed be seen in Fig. El Due to this behaviour the contribution below the lightcone 
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Wilson staggered 




40 80 40 80 



(0/T 

Figure 3: Pseudoscalar spectral functions pps{^, p)/^"^ as a function oiu /T for A''^ = 
24,32,40 and ^ = 1 (above) and vector spectral functions p)/ci;^ as a function 

of uo/T for (iV^, = (16, 1), (32, 2), (64, 4) (below). In both cases p/T = 4, m/T = 1. 



is visible in the same plot. In this case it appears that staggered fermions reproduce 
the continuum result substantially better than Wilson fermions. 

Finally, we note that the following behaviour of the kernel and spectral functions 
(p°° and pg" excluded) 

2T 

u 

uj^oo: Ph{uj,p)^uj^ A'(r, cj) ^ e""^" + e'^""^/^), (55) 

makes it difficult to study spectral functions for both small and large energies in one 
plot. This can be circumvented by instead showing the integrand at the midpoint 
r = 1/2T, i.e., 

A-(l/2r,.W(..p) = £;|^. (56) 
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Figure 4: Spectral functions p^^{u!,p)/T'^ as a function of u/T for Nr = 16 (above). 
Spectral functions Pv,a(i^, p)/?"^ sinli(a;/2T) as a function of ur/T for A^^ = 24 (be- 
low). In both cases p/T = 4, m/T = 1 and C, = 1- 



which takes a finite value for — ^ and vanishes exponentially for large u. In Fig. 
13] we show an example of this for pv,A- Since the region with large ut is exponentially 
suppressed, we note that the lattice artefacts related to the finiteness of the Brillioun 
zone discussed above give an exponentially small contribution. We conclude therefore 
that the euclidean correlator at the midpoint Gff(l/2T, p) is largely insensitive to 
these artefacts. We also point out that it follows from the analysis in Section ESI that 
the area under the curves are identical: the larger spectral weight of pA below the 
lightcone is exactly compensated by the larger spectral weight of pv above threshold. 
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5 Summary 



We have studied meson spectral functions at nonzero momentum in the infinite tem- 
perature hmit, in the continuum and on the lattice using Wilson and staggered 
fermions. We found that for large values of the energy u;, lattice spectral func- 
tions become sensitive to the effects of discretizaton and deviate from the continuum 
expectation, in agreement with the conclusions from Ref. fU]- For smaller oo, finite 
discretization affects predominantly the mismatch between the continuum and lat- 
tice lightcone, which can be substantial. In the free field limit a simple relationship 
between the euclidean correlators in different channels at the midpoint was found. 

A qualitative comparison between the results obtained with staggered and Wilson 
fermions suggests that in the low-energy region lattice artefacts are less prominent 
for the staggered formulation. The use of an anisotropic lattice, on the other hand, 
seems to be more beneficial for Wilson fermions. 
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